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ABSTRACT: Today's power system continues to expand exponentially due to extensive automated 
need for both domestic and industrial needs. The expansion of power transmission networks is limited 
due to environmental right of way and economic reasons. The ever growing demand of power system 
network has resulted in rotor angle oscillations and voltage instability. It is precisely due to the 
inability to generate the required reactive power and real power control that propounds the cause of 
voltage collapse and rotor oscillations. Therefore it is time to focus on the use of the shunt reactive 
compensation devices to alleviate this issue of voltage stability and series reactive compensation 
devices to improve the rotor angle stability. Flexible AC Transmission Systems (FACTS), which have 
recently been developed in power electronics, can make it easier to control the flow of power, boost 
power transfer capabilities, lower the cost of generating, improve security, and enhance the stability of 
power systems. The growth of large interconnected power systems demands faster response to 
maintain stability. Therefore it is also necessary to develop computationally intelligent controllers 
based FACTS devices for faster and more accurate control action. The objective of this paper is to 
analyses stability enhancements in power system networks by FACTS devices like GCSC and TSSC. In 
this work small signal stability and transient stability analysis of the above FACTS devices are 
investigated. From the analysis it is proved that TSSC is the best FACTS device for transient stability 
enhancements. 
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INTRODUCTION 


The generator rotor is provided with damper 
windings to dampen the amplitude of 
oscillations. Fast excitation systems were then 
implemented to enhance the generators’ 
reaction and maintain synchronism, although 
this decreased the oscillations’ damping. 
Therefore, it was common practice to dampen 
these oscillations using excitation systems with 
Power System Stabilizers (PSS). In order to 
increase the reliability of the power supply and 
for economic reasons, the power stations were 
heavily interconnected during the 1960s. But 
these interconnections resulted in low- 
frequency oscillations (0.1-1 Hz) in which the 
generators of one area swing against the 
generators of neighboring areas. These low 
frequency oscillations prevented further 
development of highly interconnected 
networks. Since 1960s numerous theoretical 
studies and field tests have been conducted by 
power engineers regarding the stability issues 
[1-2]. According to the nature of power system 
oscillations, controllers are developed to 
increase the damping torque. The damping of 


power oscillations is possible by using power 
system stabilizers and FACTS devices. Power 
system stabilizers are incorporated in generator 
excitation systems and are used to damp the 
generator electro-mechanical oscillations by 
taking generator speed and electrical power as 
inputs. Whereas the FACTS devices are installed 
in transmission systems away from the 
generation systems is to damp the generator 
electro-mechanical oscillations. The electric 
power generated by the generators is consumed 
by the loads at the consumer side through the 
transmission system. The transmission lines are 
operated below their thermal ratings to ensure 
reliability of electric supply to the loads. There 
is a huge increase in power demand in recent 
years due to industrialization and urbanization. 
Therefore the transmission system capacity is 
to be increased [3-5]. Due to high investment 
costs and factors like time taken for 
construction and disruption of existing power 
systems, the installation of new transmission 
lines becomes difficult. Therefore it becomes 
necessary to utilize the existing transmission 
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lines to the maximum capacity. FACTS can 
efficiently utilize the existing transmission to 
improve voltage stability of the bus and also 
damps the low frequency oscillations. The 
objective of FACTS controllers is to control the 
power flow in the transmission lines and to 
increase the transmission system reliability and 
security. The rapid development of modern 
power electronics technology has led to a 
number of FACTS controllers like GCSC and 
TSSC. Thus with marginal investment for FACTS 
device installation compared to new 
transmission line installations, it is possible to 
achieve increased power transmission within 
short period of time. Thus the FACTS devices 
installation in transmission system results in 
transmission line installation cost reduction, 
increase in transmission line power transfer 
capacity and increase in stability, security and 
reliability of power systems [6]. 


PROBLEM STATEMENT 


The rapid increase in electric power demand 
led to the burden over the electrical network 
resulting the different types of problems in the 
network. Therefore, enhancement of power 
system stability is the ultimate requirement. 
The ever growing demand of power system 
network has resulted in rotor angle oscillations 
and voltage instability. It is precisely due to the 
inability to generate the required reactive 
power and real power control that propounds 
the cause of voltage collapse and rotor 
oscillations. For improving the stability, 
different types of FACTS devices have been 
used. Various studies have been done in the 
field of power system design, modeling, 
operation, control, and application of FACTS 
devices. The detailed studies of several FACTS 
devices are strongly needed. This is because the 
operation a FACTS device affects the multiple 
conditions (contingencies, equipment failure, 
faults, sustained oscillations, etc.) occurring on a 
power system. Flexible AC Transmission 
Systems (FACTS), which can facilitate the 
control of power flow, increase the power 
transfer capability, decrease the generation 
cost, improve the security and enhance the 
stability of power systems. 


METHODOLOGY 


The methodology deals with power angle 
curve variations in the SMIB system with 
various FACTS devices for the investigation of 
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system and thus improves the dynamic stability 
of the power system. FACTS devices can 
transient stability and small signal stability. For 
analysis, a basic SMIB system is taken into 
account. The SMIB system is apply to the FACTS 
devices, such as TSSC and GCSC. Different 
FACTS devices are used to identify the systems’ 
Eigen values. Numerous inferences are seen and 
observed based on the Eigen values. The two 
scenarios mentioned above are investigated for 
TSSC and GCSC to improve power system's 
small signal stability. The conclusions reached 
showed that various FACTS devices led to 
various steady-state rotor angles and various 
rotor oscillation frequencies. If the rotor angle 
changes by just 10 degrees, the rotor 
oscillations alter until the rotor reaches its 
steady state value, and the accompanying 
changes in the There is evidence of oscillation 
frequency. In the second instance, a 0.2 pu 
change in power input is specified as a 
disturbance, leading to the observation of rotor 
oscillations and associated frequency 
oscillations. 

In both cases all FACTS devices result 
different system matrices and Eigen values as a 
result of all FACTS devices. The results received 
vary from one FACTS device to another. All of 
the cases are thoroughly discussed along with 
various stability parameters, such as the 
maximum power that can be used, the steady 
state rotor angle, the synchronizing coefficient, 
the oscillation frequency of the rotor, the 
damping ratio, the damped frequency of the 
oscillations, the response time constant, the 
system matrix, and the Eigen values. As a result, 
the analysis of tiny signal and transient stability 
of various FACTS devices integrated into the 
SMIB system is performed, and various 
conclusions are addressed based on the results. 
The SMIB system's power angle curve 
variations are tested both with and without 
FACTS devices. The small signal stability is then 
tested both the case with and without FACTS 
devices. The transient stability analysis is then 
performed both with and without FACTS 
devices. 


4. Transient Stability Analysis Using Equal 
Area Criterion 


When a machine is subjected to a significant 
disturbance, such as the sudden application of a 
heavy load, the loss of generation, the loss of a 
heavy load, or a system fault, transient stability 
studies are used to determine whether or not 
synchronism is maintained. Large disturbances 
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cause oscillations of such size that the nonlinear 
swing equation cannot be linearized and must 
instead be solved. For the purpose of predicting 
stability for significant disturbances, the Equal 
area criterion is applied. This approach relies on 
a graphical representation of the energy 
contained in the rotating mass to assess the 
machine's stability following significant 
disturbance. The following is the swing 
equation: 


H/mfo0.d26/dt2=Pm-Pe=Po (4.1) 


Where Pa is the accelerating power, H inertia 
constant, fo initial frequency, Pm mechanical 
input power, Pe electrical output power 
therefore the above equation is rewritten as 
follows 


d2 6 /dt?=2fo/H.(Pm-Pe) (4.2) 
Pm -mechanical input ,H- inertia constant, 


qô/dt change in rotor angle with respect to 
time 


The above equation after simulation results as 
follows 


dô/dt= |222, ff Pm— Peds (4.3) 
H 60 


Pe 


0 ôo ôi Ômar 
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The above equation gives the relative speed of 
the machine concerning the synchronously 
revolving reference frame for stability this 
speed must become zero at some time after 
disturbance .therefore from the above equation, 
for stability criteria the following equations are 
true 


f n Pm — Ped ô=0 (4.4.) 


A machine with initial operating angle dois 
considered for which the corresponding 
mechanical power input Pmo=Peo when there 
is a Sudden step increase in input power to 
Pm1 since Pmi>Peo the rotor starts 
accelerating and the power angle 6 increase the 
excess energy stored in the rotor during initial 
accelerating is as follows 


fÈ Pmi — Ped ô 


=area abc=area A1 (4.5) Then the energy 
stored in the rotor is given up and therefore it 
decelerated back to synchronous speed is given 
by 

fE 


sı Pm — Pe d ô =area bdc=area A2 (4.6) 


T 


Fig.1: Equal area criteria for sudden change in load 


A. Equal area criterion for sudden change in 
load and three-phase fault 


The damping contained in the machine will 
reduce these oscillations, and a new steady state 
rotor angle is 1 is reached (Fig.1). The rotor 
swings at point b between points a and d which 
is Omax. The rotor angle oscillates back and 
forth between is o and max at its natural 
frequency. This is referred to as the equal area 
criteria. Two different situations are used in the 
transient stability analysis. 1. Use in response to 


an abrupt increase in power input 2. Use with a 
three-phase fault. 


B. Application to sudden increase in power 
input 


The equal area criteria is used to determine 
the maximum additional power Pm which can 
be applied to maintain stability for a sudden 
change in power in power input the stability is 
maintained only if area A2 is less than or at least 
equal to area A1 as shown in Fig.2 the is max is 
marked for reference in Fig.2 applying the equal 
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area criteria and simplifying the equation the following expression is obtained: 


Fig.2: Equal area criterion maximum power limit 


(Smax- is 60)sin dbmax+cosémax=cosdo (4.7) stability limit can be found from 


The above nonlinear algebraic equation can be Poe Pmaxsmod (45) 


solved by iterative technique for dmax using 61=7- max (4.9) critical angle change with 
netons raphson method once max is obtained respect to change in power. 
the maximum permissible power or transient 


5. RESULTS AND DISCUSSION 


Small signal stability analysis and transient and GCSC, based on Eigen value and the system 
stability will discuss in the result part. Table 1 matrix stability of system will be verified. 
shows a small signal stability analysis of TSSC 


Table 1: Results obtained for small signal stability analysis with different FACTS devices 


1.1700pu 


1.5385pu 
0.5385rad 0.4006rad 
1.0044pu 1.4166pu 
4.3634rad 5.183 1rad 
0.2998pu 0.2525pu 
4.1636 rad 5.0152 rad 
0.6627 0.7982 
Tau 0.7643 0.7643 
th 1.3561 1.4056 
A= system 0 1.0000 0 
matrix - 37.7061 -2.6169 1.0000 
-62.3161 - 
2.6169 
K=eign value -1.3085+5.9995i- | -1.3085 + 
1.3085 - 5.9995i 7.78491 - 
1.3085 - 
7.78491 
The eigen values of A has a negative real part same as the negative real part of the Eigen 
for both cases and therefore lie on the left half values are same whereas there is a slight 
of the complex plane, hence the system is stable. increase in the frequency of oscillations for 


In both, the cases the damping remains the TCPAR-based SMIB. Thus the TCPAR improves 
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the system stability by increasing the reducing the power angle, 5. The various values 
synchronizing power constant and thereby obtained are shown in Table 1. 


A Case study for Transient stability analysis 
5.1.1 Case 01 Sudden Increase in power input with FACTS device 
Sudden Increase in power input with GCSC (Pmo =0.6) 


a 


Power angle curve GCSC 


aT T : 
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Fig.3 MATLAB plot of equal area criterion with GCSC for a sudden change in load 


With GCSC connected to it, the SMIB system greater than 126.460° degrees. For the rotor to 
under consideration is impacted by a sudden lose synchronism and which leads to an 
shift in power input. According to the equal area unstable condition. Thus it is inferred that to 
criterion curves plotted in MATLAB as shown in maintain stability the sudden increase in power 
Fig.3, the mechanical power input should not be input should be less than the maximum limit 
greater than 1.977 pu, which is referred to as otherwise which results in a loss of 
the total power of critical stability, for a synchronism with the infinite bus. And thus an 
machine operating with an initial power input additional sudden power input of 0.293 pu can 
of Pm = 0.6 pu and an initial operating angle of be given for the SMIB system with GCSC. There 
15.466 degrees. The sudden maximum is a increase in the maximum angle of swing 
additional power that can be applied is 1.210 from 125.840° to 126.460. Thus it is clear that 
pu. A new operating angle of 53.5400 is then GCSC improves the transient stability of the 
reached. The rotor will lose synchronism and system by increasing the transient stability 
become unstable if the maximum swing angle is margin. 


Case 2 Sudden Increase in power input with TSSC (Pmo = 0.6): 


Power angle curve TSSC Facts device 


Electrical power 
© 
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0.4F 


0 1 1 1 L 1 L 1 1 
0 20 40 60 80 100 120 140 160 180 
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Fig.4: MATLAB plot of equal area criterion with TSSC for sudden change in load 


The SMIB system considered is subjected to curves plotted in MATLAB as shown in Fig.4, it 
the sudden change in power input with TSSC is observed that for a machine operating with 
connected to it. From the equal area criterion an initial power input of Pm = 0.6 pu, initial 
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operating angle = 12.840 degrees, the sudden 
maximum additional power that can be applied 
is 1.537pu that is the mechanical power input 
should not exceed 1.6 pu, which is called as the 
total power of critical stability. Then it settles at 
a new operating angle ô= 52.317 The rotor will 
lose synchronism and become unstable if the 
maximum swing angle is greater than 127.68° 
Thus, it can be concluded that to retain stability, 
a rapid increase in power input must be greater 
than the maximum allowable value; otherwise, 
synchronism with the infinite bus is lost. The 
Table 2 compares the parameters of the SMIB 
transient stability system with and without 
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TSSC, and it is clear from this comparison that 
the SMIB with TSSC increases the transient 
stability margin by decreasing the initial power 
angle by margin by increase the initial power 
angle by 1.537. It is also observed that sudden 
additional power with TSSC is increase from 
1.084 pu to 1.537 p.u. and thus an additional 
sudden power input of 0.235 pu can be given for 
the SMIB system with TSSC There is a increase 
in maximum angle of swing from 125.840° to 
127.68 deg. Thus it is clear that TSSC improves 
the transient stability of the system by 
increasing the transient stability margin. 


Table 2: Results obtained for Transient stability analysis (Sudden Change in load) with different FACTS devices 


| Initial power | 0.600pu | 0.60pu | 0.600pu 
oe 16.7919 15.466° 12.840° 
Initial power angle 
Sudden additional power 1.084pu 1.210pu 1.537pu 
Total power for critical stability | 1.684pu 1.810pu 2.137pu 
Maximum angle swing 125.8409 | 126.460° | 127.683° 
New operating angle 54.160° 53.5400 52.3179 


5.3 Three Phase fault with FACTS devices 
Sending end and Midpoint fault with GCSC 


GCSC fault at sending end 


Electrical power in pu 
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Fig.5: MATLAB plot of equal area criterion with GCSC for 3 phase fault at sending end 


The transmitting end three-phase fault near 
the generator is present in the SMIB system 
taken into account by the GCSC. For a machine 
working at Pm = 0.8 pu, with an initial operating 
angle of 6.807 degrees, it is noted using the 
equal area criterion plotted in MATLAB as 
shown in Fig.5 that the fault should be cleared 


by the time shorter than the critical clearing 
time, tc = 0.219 sec. If a fault is fixed after a 
certain amount of time, the system loses 
synchronization with the infinite bus and 
becomes unstable. The maximum angle of swing 
at tc is the angle above which the rotor losses its 
synchronism which is equal to 154.612°. If the 
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fault is not cleared before critical clearing angle 
75.93 deg then the accelerating energy area is 
increased such that decelerating area cannot be 


fault at midpoint(gcsc) 
T 
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greater than the accelerating energy area and 
therefore the rotor loses its synchronism with 
infinite bus and the machine becomes unstable. 


a 
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Fig.6: MATLAB plot of equal area criterion with GCSC for 3 phase fault at mid-point 
of transmission line 


The SMIB system considered with GCSC is 
subjected to midpoint three phase faults 
between the generator and the infinite bus. 
Using the equal area criterion plotted in 
MATLAB as shown in Fig.6 it is observed that 
for a machine operating at Pm = 0.8 pu, with 
initial operating angle 6.807 degrees, the fault 
should be cleared within the rotor reaches the 
critical clearing angle 120.465 deg. 


The system loses synchronism with the 
infinite bus and becomes unstable if the 


Sending end and Midpoint fault with TSSC 


25 


problem is fixed after the rotor has beyond the 
critical clearing angle. The critical clearing 
angle, which is equivalent to 120.465, has a 
maximum swing angle over which the rotor 
loses synchronism. The rotor loses synchronism 
with the infinite bus and the machine becomes 
unstable if the defect is not fixed before the 
critical clearing angle of 75.93, which causes the 
accelerating energy area to increase such that 
the decelerating area cannot be bigger than the 
accelerating energy area. 


Electrical Power 


0 20 40 60 80 


Rotor angle in deg 


Fig.7: MATLAB plot of equal area criterion with TSSC for 3 phase fault at sending end 


The SMIB system considered with TSSC is 
subjected to sending end three-phase fault near 
the generator. Using equal area criterion plotted 
in MATLAB as shown in Fig.7,it is observed that 
for a machine operating at Pm = 0.8 pu, with an 
initial operating angle 17.235 degrees, the fault 


should be cleared by the time less than critical 
clearing time, tc = 2.298 sec. If the fault is 
cleared after critical clearing time then the 
system losses its synchronism with the infinite 
bus and becomes unstable. The maximum angle 
of swing at tc is the angle above which the rotor 
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losses its synchronism which is equal to 
162.765°. If the fault is not cleared before 
critical clearing angle 101.684° then the 
accelerating energy area is increased such that 
decelerating area cannot be greater than the 
accelerating energy area and therefore the rotor 
loses its synchronism with the infinite bus and 
the machine becomes unstable. 


Equal area criterion-Mid Point 3 phase fault-- 
with TSSC connected 


Three phase midpoint faults between the 
generator and the infinite bus are present in the 
SMIB system taken into account using GCSC. For 
a machine working at Pm=0.8 pu, with a 
starting operating angle of 17.235 degrees, it is 
noted using the equal area criterion plotted in 
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MATLAB as shown in Fig.8 that the fault should 
be cleared before the rotor reaches the critical 
clearing angle of 141.13 deg. If the fault is 
cleared after the rotor exceeding the critical 
clearing angle then the system losses its 
synchronism with the infinite bus and becomes 
unstable. The maximum angle of swing at 
critical clearing angle is the angle above which 
rotor losses its synchronism which is equal to 
162.71°. If the fault is not cleared before critical 
clearing angle 146.21°, then the accelerating 
energy area is increased such that decelerating 
area cannot be greater than the accelerating 
energy area and therefore the rotor loses its 
synchronism with infinite bus and the machine 
becomes unstable. 


= 
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Electrical power 
[e] 
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ie) 20 40 60 
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Fig.8: MATLAB plot of equal area criterion with TSSC for 3 phase fault 
at midpoint of transmission line 
Table 3 Results obtained for Transient stability analysis (Three phase fault) with 
different FACTS devices 


NO , 
FACTS Sending end fault | 22.665 | 157.345 90.99 0.397 
device Midpoint fault 22.665 | 155.492 107.86 

Sending end fault | 6.807 | 154.612 75.938 0.419 
GCSC 

Midpoint fault 6.807 | 153.612 120.465 

Sending end fault | 17.235 | 162.765 101.684 | 0.498 
TSSC 

Midpoint fault 17.235 | 162.765 141.381 


5.4 Critical Clearing Time for Mid-Point Fault 
Using Swing Curve without FACTS 


For mid-point fault the critical clearing time 
could not be found out using the equal area 
criterion, because of the nonlinear swing 


equation. To determine the crucial clearance 
time for mid-point fault, the nonlinear swing 
equation is therefore solved using the Runge- 
Kutta method. Plots of the swing curve are 
shown in Fig.9 for both the SMIB system with 
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and without the FACTS device. The input is the 
fault clearance time. The power angle in the 
swing curve increases without bounds and the 
system becomes unstable if the fault clearing 
time is longer than the critical clearance time. In 
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contrast, when the power angle 6 returns to 
the swing curve following the greatest swing, 
the curve returns after the maximum swing 
then the oscillations will subside and the system 
returns to the stable state. 


One-machine system swing curve. Fault cleared at 0.3s and 0.5s 
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Fig.9: Swing curve without FACTS 


The SMIB system considered is subjected to a 
midpoint fault between the generator and the 
infinite bus. The power angle, 5 versus time 
called the swing curve is plotted using MATLAB 
as shown in Fig.9. The critical clearing angle 
was found to be equal to 107.869° for which the 
corresponding critical clearing time is 0.44 
seconds. The graph shows that when the fault 
clearance time, tf is less than the critical 
clearing time, tc , the system returns back to a 
stable state and synchronizes with the infinite 
bus. Whereas when the fault clearing time tf is 
greater than the critical clearing time, the 
power angle, 5 keeps on increasing beyond 
1800 , indicating that the system stability is lost 
with the infinite bus as it loses its synchronism. 


Critical clearing time for mid-point fault 
using swing curve with FACTS 


The SMIB system considered with GCSC is 
subjected to a midpoint fault between the 
generator and the infinite bus. The power angle, 
6 versus time called swing curve is plotted using 
Matla shown in Fig.10. It was found that the 
critical clearing angle is equal to 75.9380, and 
the critical clearing duration is 0.419seconds. 
The graph demonstrates that the system 
returns to a stable condition and synchronizes 
with the infinite bus when the fault clearing 
time, tf, is smaller than the critical clearing time, 
tc. 


One-machine system swing curve. Fault cleared at different time 
T T p T 7 4 


5007 


> 

[=] 

[=] 
T 


Delta, degree 
wo 
Q 
O 


N 

S 

S 
T 


100 F 


—t=0.4sec 
—t=0.5sec 
—t=0.6sec| | 
= = t=0.8sec 


0.6 


0.8 1 


t, sec 


Fig.10: Swing curve with GCSC 
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Whereas when the fault clearing time tf is 
greater than the critical clearing time, the 


power angle, 5 keeps on increasing beyond 
400 T T 
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180°, indicating that the system stability is lost 
with the infinite bus as it loses its synchronism. 
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Fig.11: Swing curve with TSSC 


The SMIB system considered is subjected to a 
midpoint fault between the generator and the 
infinite bus. The power angle, 5 versus time 
called swing curve is plotted using MATLAB as 
shown in Fig.11. The critical clearing angle was 
found to be equal to 101.684 for which the 
corresponding critical clearing time is 0.498 
seconds. The graph shows that when the fault 


clearance time, tf is less than the critical 
clearing time, tc , the system the system returns 
back to a stable state and synchronizes with the 
infinite bus. Whereas when the fault clearing 
time tf is greater than the critical clearing time, 
the power angle, 6 keeps on increasing beyond 
180°, indicating that the system stability is lost 
with the infinite bus as it loses its synchronism. 


Table 4 Comparison of all FACTS devices for - Swing curve and critical clearing time 


107.869 0.392 


01 No facts device 
02 GCSC 75.938 0.419 
03 TSSC 101.684 0.496 


From the Table 4 the inference drawn is the 
facts device which having higher critical 
clearing time which gives an better stability 
limit because system returns back to a stable 
state and synchronizes with the infinite bus. 
Hence TSSC is best FACTS device compared 
with GCSC. 


6. CONCLUSION 


A comparison of various FACTS devices, such 
as the GCSC and TSSC, is offered for improving 
small signal and transient stability. By 
computing the system matrix and its 
corresponding Eigen values, small signal 
stability is carried out using Eigen value 
analysis, and several implications for improving 
small signal stability are presented for various 
FACTS controllers. When several FACTS 
controllers are coupled to the SMIB system, the 
transient stability experiments are conducted 
using the equal area criterion. For the SMIB 
system controlled by FACTS, parameters like 
the critical clearing angle and critical clearing 


time are established. It is evident from the 
analysis that impedance base FACTS controllers 
are offering more stability. In contrast to 
alternative FACTS controllers, it is also obvious 
that TSSC offers greater stability. Small-signal 
stability studies of transient stability in small 
and transient stability studies the computation 
and analysis with respect to sending the end 
bus can only be made. That is the dynamic and 
fault analysis on sending end bus of the power 
system can only be performed. Whereas, the 
load side transmission system parameters 
cannot be computed and analyzed. 
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